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Abstract 

We introduce the quantum complexity class FQMA. This class describes the com- 
plexity of generating a quantum state that serves as a witness for a given QMA problem. 
In a certain sense, FQMA is the quantum analogue of FNP (function problems asso- 
ciated with NP). The latter describes the complexity of finding a succinct proof for a 
NP decision problem. Whereas all FNP problems can be reduced to NP, there is no 
obvious reduction of FQMA to QMA since the solution of FQMA is a quantum state 
and the solution of QMA the answer "yes" or "no" . 

We consider quantum state generators that get classical descriptions of 3-local 
Hamiltonians on n qubits as input and prepare low energy states for these systems 
as output. We show that such state generators can be used to prepare witnesses for 
QMA problems. Hence low energy state preparation is FQMA-complete. Our proofs 
are extensions of the proofs by Kitaev et al. and Kempe and Regev for the QMA- 
completeness of fc-local Hamiltonian problems. We show that FQMA can be solved by 
preparing thermal equilibrium states with an appropriate temperature decreasing as 
the reciprocal of a polynomial in n. 



1 Introduction 



In analogy to classical computer science one may describe quantum complexity theory 
as the field dealing with the complexity of solving computational problems when a 
quantum computer is available. However, in contrast to classical algorithms, the task 
of a "quantum algorithm" is not necessarily to solve a computational problem. It 
may, for instance, also be the task to prepare a certain multi-particle quantum state 
which may be a resource for various applications 1 . Since the preparation of generic 
entangled states in multipartite systems is non-trivial, it makes sense to consider the 
complexity of state preparation procedures |2j. 
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Remarkably, the question of complexity of state preparation comes up straightfor- 
wardly when some computational problems are generalized to the quantum domain. 
This shows the following example: Consider a function 

/: {0,1}* -{0,1} 

that is given by a boolean expression of polynomial size in n. The task is to decide 
whether there exists an input x such that f(x) = 1. This problem is called SAT 
(Satisfiability). It is NP-complete j^j. The problem to find the input x whenever it 
exists is called FSAT. In general one defines FNP as the problem class that consists of 
the problems to find proofs to NP decision problem. It is easy to see that FSAT can 
be reduced to SAT (with polynomial overhead) as follows: Define the functions f° and 
f 1 that are obtained if the first input bit is set to or 1, respectively and ask whether 
there exists an x' of size n — 1 such that / has the output 1 or that f 1 has (see [3]). 
This determines the first bit of the solution. In the same way the other bits can be 
determined with a polynomial number of calls of the decision problem 1 . 

The quantum analogue of this setting is as follows. Consider a unitary quantum 
circuit U of polynomial size (i.e. consisting of a polynomial number of two-qubit gates) 
acting on 

( C 2)®(n+m) where 

we have n input and m ancilla qubits. The ancilla register 
is initialized to the state |0 . . . 0) in order to mimic irreversible mappings such as, for 
instance, the classical logical functions AND, NAND, NOR, OR [I]. According to [3] 
the quantum analogue to SAT is as follows. Given a quantum circuit U that acts on 
a tensor product of an input state \tp) and the state |0...0) of an initialized ancilla 
register. We obtain 

|0) :=17(|^>®|0...0» 

and measure the right most qubit. Then the quantum analogue to "expression is not 
satisfiable" is that for all states the measurement result is with high probability 
"no". The analogue to "satisfiable" is that there is a state such that the result is 
likely to be "yes" . This class of problems is referred to as BQNP [5j or QMA jSJ [7] 
and each state which leads to "yes" with high probability is called "witness" of the 
QMA-problem. 

Note that the problem class QMA is only to decide whether there exists a witness 
and not to prepare it if it exists. It should be noted that we see no obvious reduction 
of the preparation problem to the decision problem QMA. Consider the case that there 
is an entangled state \ip) which is accepted by U with high probability. This does not 
imply that there is any state of the first qubit of the input register such that 
IV'i) ® W) i s accepted for an appropriate state Therefore, already the first step 
of the reduction above fails. This may be a substantial difference between quantum 
and classical: The class FNP given by the set of problems to find a witness reduces 
to NP. The problem to prepare a "witness" of a QMA-problem shall be called FQMA 
here. Note that it is not even clear whether there exists a polynomial description of the 
witness state \tp). Therefore FQMA is, in the definition proposed here, not to describe 
the state that is accepted, we demand rather to prepare it. In other words, the answer 
of a QMA problem is classical whereas the answer of FQMA is quantum. 

1 Note that there are FNP-problems which are not canonically associated with NP-problems. They can 
be reduced to appropriate NP-problems nevertheless [3]. 
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In this paper we shall provide an example of an FQMA-complete problem which 
has a rather intuitive physical meaning. It is the task to prepare low-energy states in 
many-particle systems. This problem is strongly related to the invention of efficient 
cooling mechanisms as we will discuss later. The relation between efficient cooling of 
many-particle systems with non-trivial interactions and solving hard computational 
problems is well-known (see e.g. ■SJ. For Ising-type interactions one can show that 
finding ground states requires to solve NP-complete problems However, these are 
examples where the solution is only a classical spin configuration which is easy to 
describe. A natural question is the following. Consider Hamiltonians for which it is 
not clear whether a short classical description exists, is it a hard quantum control 
problem to prepare low energy states? As usually in complexity theory, we can only 
show that it is as hard as another problem, namely solving FQMA in general. 

2 Introducing the complexity class FQMA 

In order to motivate our definition, we recall the definition of NP and FNP problems. 
A language, i.e., subset of the set of strings is in NP if there exists a polynomial- 
time decidable, polynomially balanced relation Rl such that there is a string y with 
Rl(x,u) if and only if x £ L. The function problem associated with L, written FL is 
the following computational problem: 

Given x, return a string y such that Ri,(x,y) provided that such a string exists; if 
no such string exists, return "no". Less formally it is the problem to give a proof that 
the answer is "yes" whenever it is true. 

The complexity class QMA consists of the problems to decide whether a given string 
is in a certain language in QMA. The set of QMA languages is defined following 7 . 
First we introduce some notations that will be used trough the paper. We denote the 
Hilbert space of a qubit by B := C 2 . Let x S {0, 1}* be an arbitrary binary string. We 
denote the length of x by \x\. For any Hilbert space S(TC) we denote the set of density 
matrices acting on H by S(7i). 

Definition 1 (QMA) 

Fix e = e(\x\) such that 2~^ x ^ < e < 1/3. Then a language L is in QMA if for every 
classical input x 6 {0,1}* one can efficiently generate (by classical precomputation) a 
quantum circuit U x ("verifier") consisting of at most p(\x\) elementary gates for an 
appropriate polynomial p such that U x acts on the Hilbert space 

where n x and m x grow at most polynomially in \x\. The first part is the input register 
and the second is the ancilla register. Furthermore U x has the property that 

1. For all x £ L there exists a quantum state p that is accepted by the circuit with 
high probability, i.e., 

3pe S{B nx ), tr(U x (p ® |0...0)(0...0|)C/| Pi) > 1-e, 

where Pi is the projection corresponding to the measurement "Is the first qubit in 
state 1 ?". 
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2. For all x L all quantum states are rejected with high probability, i.e., 

VpG S{B n *), tr(U x (p ® |0 . . .0)(0 . . .0|)E4 -Pi) < e. 

Note that our "proofs" are mixed states in contrast to the definitions in jSJ [7j . Due 
to linearity arguments this modification does not change the language L. Note fur- 
thermore that it is always possible to construct a verifier for the same language with 
e' arbitrarily close to 0. This "amplification of probabilities" is described in in full 
detail. Since this amplification procedure has to be modified for FQMA-problems we 
briefly sketch the idea. Use k = poly(\x\) copies of the circuit U x . The decision whether 
the state p' on (B nx <&B mx )® k is accepted is based on a kind of majority function consid- 
ering all outcomes of measurements on the first qubits of all the k copies of the original 
register. We say that the total circuit accepts p' if at least k/2 measurements return 
"yes" and rejects otherwise. Actually, this decision procedure does not fit strictly in 
the setting above since it relies on k measurements. This inconsistency can easily be 
removed by realizing the majority vote by a "small" additional quantum circuit. As 
shown in [2] there exists a joint state on the k copies of the original register that is 
accepted with high probability if and only if there exists a state on one copy that is 
likely to be accepted. 

Now we define FQMA. 

Definition 2 (FQMA) 

A channel with classical input x and quantum output p x is in FQMA if there is a 
language L in QMA with a verifier U x and e as in Definition^ such that 

tr(U x (p x ®\0...0)(0...0\)U x iP 1 ) >l-e-5, 

whenever x G L, where 5 = l/r(|x|) for an arbitrary polynomial r. For x G" L the 
output is allowed to be arbitrary. 

3 Low energy state preparation 
is FQMA-complete 

The determination of eigenvalues of non-trivial many-particle interaction Hamiltonians 
is known as a computationally hard problem. This is even the case if the Hamiltonians 
are restricted to A:-local ones. Here we call a Hamiltonian on n qubits /c-local if it is 
a sum of operators that act only on k qubits. Already for 2-local interactions one has 
NP-completeness ,9.. For 3-local terms it has recently been shown by Kempe and Regev 
[7] (as extension of theorems by Kitaev et al. jSJ) that it is QMA-complete to decide 
whether the Hamiltonian has an eigenvalue less than a given value. More explicitly, 
they constructed a 3-local Hamiltonian associated with a circuit U (consisting of L 
two-qubit gates) which has an eigenvalue less than 

e 

L + l 



4 



if case 1 in Definition ^ holds. Conversely, all eigenvalues are at least 



c 



L 3 



if case 2 holds, where c is a constant. Moreover, Kitaev et al. have shown that 
the problem to decide whether a fe-local Hamiltonian has an eigenvalue below a or all 
eigenvalues are above b is in QMA whenever the gap b — a is not stronger decreasing 
than polynomially. 

In the following we will show that the arguments in [7] and [3] can be extended to 
prove that preparing low energy states is FQMA-complete. Explicitly, we define this 
problem class as follows. 

Definition 3 (fc-local Low Energy State Preparation) 

Let (H n ) be a sequence of k-local interaction Hamiltonians acting on B n . The terms 
in the sum of each H n are assumed to be positive and of norm not greater than 1. 
Furthermore, they are specified with p(n) bits where p is a polynomial. Let (a n ), (b n ) 
and (d n ) be sequences of numbers such that b n > d n > a n and the gaps between b n and 
d n and between d n and a n decrease only polynomially in n. It is promised that either 
all eigenvalues of H n are greater than b n or there exists an eigenvalue smaller than a n . 

Then Low Energy State Preparation is the problem to prepare density matrices p n 
on n qubits (after H n has been specified by a classical input) satisfying the low energy 
condition 



for all n with H n having an eigenvalue less or equal than a n . 
Lemma 1 The problem k-local Low Energy Preparation is in FQMA. 

Proof: As Kitaev et al. |2j have shown there exists a language in QMA such that 
x £ L if there exist an eigenvalue smaller or equal than a n and x £ L if all eigenvalues 
are greater or equal than b n . We have to rephrase the verifier in order to show that a 
state generator preparing witness states is also able to prepare states with mean energy 
less or equal than d n . 

For fixed n we drop the index n and write H as sum of its /c-qubit terms 



Since < Hj < 1 for all terms, the operators Hj, 1—Hj define a positive operator valued 
measure. This POVM can be implemented with a constant number of quantum gates, 
ancillas and one-qubit measurements. Let the operator Hj correspond to the result 
and 1 — Hj to the result 1. Use a random generator that has one number 1, ... ,r 
as outcome with equal probability. If the outcome is j perform the measurement 
Hj , 1 — Hj. The random generator can be implemented by an appropriate quantum 
circuit and the measurements can be replaced by controlled bit-flips on an additional 
ancilla which shows the measurement result. This defines a quantum circuit U and 



tr(p n H n ) < d, 



it 
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the latter ancilla qubit is the qubit that is measured. Thus we have the setting of 
Definition [2 Then the probability that a state p is rejected is given by 



Assume there exists a state with energy at most a n . Then it is accepted by the circuit 
with probability 1 — a n /r. A channel which is a "universal FQMA solver" is able to 
prepare a state which is accepted with probability at least 1 — d n /r, i.e., with mean 
energy at most d n . This completes the proof. □ 

As we will prove in the following, we have also the other direction: 

Theorem 1 (FQMA-completeness) 

The Problem 3-local Low Energy Preparation is FQMA- complete, i.e., each quantum 
state generator that solves the 3-local Low Energy Preparation Problem can be used to 
generate witnesses for QMA problems with arbitrary precision. Whenever there exists 
a state that is accepted by a given circuit with probability 1 — e the "Low Energy State 
Generator" can be used to prepare a state that is accepted with probability 1 — e — 5 
where 5 can be made arbitrarily small. The complexity of the procedure is 0(1/8 ). 

Given a quantum circuit U (note that we shall drop the subscript x in the following) , 
the task is to construct a Hamiltonian such that its low energy states can efficiently be 
transformed into a state that is accepted by U with high probability. 

The constructions of Kitaev et al. jS] and Kempe and Regev |7J show a correspon- 
dence between states that are accepted by a given circuit and low energy states of the 
corresponding Hamiltonian. One may therefore expect that our Theorem is already 
given by a straightforward reinterpretation of their results. However, the main problem 
is that such a straightforward approach allows only to prepare states that are accepted 
with considerably lower probability than 1 — e. The QMA problem is to decide whether 
there exists a state that is accepted with probability at least 1 — e or whether all states 
are accepted with probability at most e. The gap between e and 1 — e is essential for the 
proofs in d • We cannot make use of such a gap. Whenever there exist states that 
are accepted with probability 1 — e and others with probability e there are, by linearity, 
always states that are accepted with probability p for every value between 1 — e and e. 
In the straightforward extension of [3] and [7j, it is not clear how to avoid that some 
of these states are obtained. Therefore we demand only to construct states which are 
accepted with probability 1 — e — S. As we will see below, the allowed inaccuracy 5 is 
essential for our construction. 

The main idea is to construct a circuit U which is a modification of the amplification 
explained in Section |2j For the "meta-circuit" U exist states that are accepted with 
probability 1— e where e is considerably smaller than e. Now we can use the construction 
of [7] and obtain a Hamiltonian. We show that its low energy states can be transformed 
simply to states are accepted by U with high probability. Although the probability of 
acceptance may in general be considerably smaller than the possible value 1 — e we can 
use them to obtain states which are accepted by the original circuit U with probability 
almost 1 — e. 
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To motivate our construction of U one should note that Kitaev's probability ampli- 
fication does not provide states that are accepted by the original circuit U with high 
probability. Assume, for instance, that we use k = 31 copies of U for large I G N. Let 
the input of U be the state 

p ml ®^\ (1) 

where p is a state that is always accepted by U and 7 is a state that is never accepted. 
The joint state (0) is very likely to be accepted by the majority vote. A straightforward 
method to obtain an input state for the original circuit would be to choose randomly 
the restriction of the joint state to one of the k input registers. We would get the state 

2 1 

This state is rejected with probability 1/3 independent of k. 

Therefore we modify the "majority vote" in such a way that we accept the joint state 
p only if more than fe(l — e — 1/ x^k) of the original circuits U accept their corresponding 
parts of p. By choosing randomly one of the k reduced states one obtains a mixed state 
that is accepted by U with conditional probability at least 

1-e-l/Vk, 

given the event that the meta circuit has accepted p. 

The following lemma shows quantitatively that there exists a state that is accepted 
by the meta circuit with high probability: 

Lemma 2 Given the modified majority vote which accepts input states of U x whenever 
more than k(l — e — 1/yk) copies ofU x have accepted. Let p be a state that is accepted 
by U x with probability at least 1 — e. Then a k-fold copy of p is rejected with probability 

0(2-^/ ln2 ), 

i.e. the probability can be amplified sub-exponentially. 

Proof: Assume the worst case that U accepts p exactly with probability 1 — e. 
In order to have a lower bound on the probability that p® k is accepted we have to 
estimate the tail of a binomial distribution according to a Bernoulli experiment with 
probability e: 

P{p® k rejected ) = £ - e) k ^ =: , 

j<i j<i 

where 

l:= \k(l-e-l/\fk)]. 

For I > k/2 we have certainly Ylj<i^k(j) < The binomial coefficient b/-(l) can 

be estimated using relative entropies [TO] : 



b k (l) = 0(2- kD (W) 
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with 



D(l/k) := (l/k) log 2 e - (l/k) log 2 (l/k) + (1 - l/k) log 2 (l - e) - (1 - l/k) log 2 (l - l/fc) 

Note that the term D(l/k) is the Kullback-Leibler relative entropy which measures 
the distance between the probability measure on two points given by (e, 1 — e) and the 
(formal) probability measure defined by the observed relative frequencies (l/k, 1 — l/k). 
Hence p® k is rejected with probability 

We use Taylor expansion around the point 1 — e up to the first order with an appropriate 
remainder term and obtain 

D (l- e - X /¥k) = D(l - e) + D'(l -e)l/<fk + L>"(0)_L_ , 

for an appropriate value G in the interval (1 — e, 1 — e + 1/ v^)- Relative entropy is 
always strictly positive for two different probability measures. Therefore the function 
D has its unique minimum at 1 — e and the first and second term in the expansion 
vanish. The term D"(Q) is at least 2/ In 2. We have 



D(l -e-l/Vk) > 



\n. 2y/k ' 
We conclude 

P(p® k rejected ) = 0(2~^' ln2 ) . 

□ 

Now we are able to prove Theorem ^ The first part of the proof rephrases the 
constructions and notations of and 0. Proof: Let U = UlUl-i ■ ■ ■ U\ consist 
of L two-qubit gates. The circuit U acts on B l = B n B m qubits, where n and m 
are the size of the input and the ancilla register, respectively. Furthermore we need 
an additional register called "clock". It consists of L qubits. The total Hamiltonian 
consists of four parts: 

" • — Hin ~\~ Hout ~\~ Hprop clock ■ 

The low energy states of H represent in some sense the whole history of the quantum 
circuit U. The correlations between the clock register and the Hilbert space the circuit 
acts on contain the information at which time step which gate has been applied. This 
can be achieved as follows. 

H c lock ■= L 12 |01)jj(01|jj . 

l<i<j<L 

The subscripts i,j indicate the considered qubits. H c i oc k acts only on the L clock 
qubits. It penalizes all states in the clock register that are not of the form 

1 11... 1^0) , 



called "unary representation" of the numbers 1,...,L. States of this form are denoted 
by \t). They are the allowed states of the clock. 
Consider i?j n . It is defined as follows. 



N 



H in := |l}i<l|i®|0)i(0|i 



i=m+l 

The first component of the tensor product acts on the ancilla register and the second 
on the first qubit of the clock. It penalizes all states where the ancilla register is not 
initialized whenever the clock is in its starting position. 

Consider H prop . It ensures that the correlations between the clock and the remain- 
ing registers are according to the history of the quantum circuit. 

L 



t=l 

with 

H P ro P ,t := \{l ® |10) M+1 (10| M+1 + 1 ® |10) t _i >t <10|t- 1)( -U t ® \l) t (0\t + Uj\0)t(l\ t ) 
for 2 < t < L - 1 and 

H P rop,i := \{l ® |10)i, 2 (10|i, 2 + 1 ® |0)i(0|i - U X (8) |l)i(0|i + E/}|0)i(l|i) 

Hpro P ,L := ^(1 ® + 1 8> |10) L _i, L (10| L _i, L -U L ® |1)l(0|l + ul\Q) L (l\ L ) . 

Consider H ou t- It penalizes all states where the output is not in the state |0) 
whenever the clock is in its end position. 

Hout ■= |0>i(0|i<g> |1) L (1| L . 
We also need unitary transformation W defined in [S] 

L 

W:=J2 U tUt-i---U 1 ® |t)(t|. 
t=o 

for transforming the low energy states. This transformation can obviously be imple- 
mented efficiently since it can be written as product of transformations 

W t :=U t ®\l) t {l\t, 

where each Wt is a 3-qubit gate. 

Having introduced the necessary notation we come to the second part of the proof. 
Let p be a state such that tr(pH) < a and a be the restriction of the state W^pW to 
the input register. We show that a is accepted by the circuit U with probability at 
least 

I - a 1 ' 4 poly(L) . 
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Here and in the rest of the proof we write poly(L) for every expression f(L) which is 
smaller than an appropriate polynomial in the size L of the circuit independent of its 
specific structure and independent of a. 

Without loss of generality we may assume p to be pure due to linearity of mean 
energy. Let W^pW be the state \4>}(4>\. We decompose \(f>) as follows 

\4>) = d \a ) <8> |-0o> + ^2dj\aj) <g> \tpj) , 

where is the j-th eigenvector of W^(H prop + H c i oc k)W (starting with 0). The only 
eigenvector corresponding to the eigenvalue is 



The other eigenvalues are bounded from below by the non-zero eigenvalues of H prop . 
They are given by = 1— cos qu with q^ := 7r/c/(L+l). Asymptotically, the eigenvalues 
are larger than c/L 2 for an appropriate positive constant c. We have 

a > tr{W ] P WW ] HW) = (<f>\W* HW\<f>) > (<p\W^ H prop W\<p) > ^ \dj\ 2 qj . 
We conclude 

^2 \dj\ 2 < a poly (L) . 
Now we neglect the terms for j > 0. The remaining state is 

It) : = l«o) <8) \tpo) ■ 

The norm distance between \<f>) and I7) can now be estimated by elementary geometry. 
Let \x) and \y) be two vectors such that \z) = \x) + \y) is a unit vector. Then we always 
have 

llk>- ip|||x>||<v^|||y>||, 

since the left hand side is the length of a hypotenuse of an rectangular triangle with 
legs of length |||y)|| and 1 — |||x)|| < |||y)||. 

Hence the norm distance between I7) and the true state vector \cj>) is less than 

\J 2 Ej>i \ d j\ 2 = Vapoly(L) 
We decompose (7) as 

|ao)®|^)=co|?to)®|0...0)®hfo)+ Yl c b\Vb) ® \b) ® \i>o) , 

l<b<2 m 

where \b) denote the basis states of the ancilla register and 1%) are arbitrary state 
vectors in the input register. 

In the following we will use the following elementary argument several times: If the 
norm distance between two unit vectors \x) and \y) is A < 1 one has 

(x\A\x) - (y\A\y)\ < \\A\\3A, 
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where A is an arbitrary matrix with operator norm \\A\\. Using this estimation we have 
a + ^poly{L)\\H\\ > {-r\W^HW\i) > (-/\W^ H m W\ 7 ) = j^— ^ \ c b\ 2 • 

l<b<2 m 

Since \\H\\ < poly(L) we find 

^2 \cb\ 2 < Vapoly(L) . 

l<b<2 m 

We neglect the part of the state I7) with q, for b > 1. The remaining state has the 
form 

Ix) := |»fc}®|0...0)® 
Its norm distance from I7) is a 1//4 pc%(L). Its mean energy is at least 

P(\X) rejected ) 



( X \W^H out W\x) 



L + l 



The norm distance between \<j)) and |x) is a 1 / 4 poly(L). By the geometry argument 
above and ||-ff OU i|| = 1 we have 

a > WWlH^WW > P( ' X) /f f 6 ' } - a^poly(L) 

Li -\- 1 

This implies 

P{\x) rejected ) < a 1 / 4 poly(L) . 

This shows that for each state \(f>) with energy a the corresponding state \x) is re- 
jected with probability a 1//4 poZy(L). Since the norm distance between \x) and \<p) is 
a^polyiL) we conclude that each state with energy smaller than a is rejected with 
probability at most 

e> :=a 1/4 po/y(L). (2) 

The Hamiltonian H corresponding to the meta circuit U as explained in the first part 
of the proof has an eigenvalue of the order 

( 2 -Vfc/ln2-j _ 

Due to [Sj and [7] this is guaranteed by the fact that U has an input that is rejected 
with this small probability. Hence we can use the Low Energy State Generator for 
preparing a state p with energy d < l/s(L) for any polynomial s. 

Let p be the restriction of p to a randomly chosen sub-register. The conditional 
probability that it is accepted by U given the event that p has been accepted by the 
meta circuit U is at least 

1 

1 - e- -77= 

as already explained above. The meta circuit accepts with probability 1 — d l / A poly{L). 
The unconditional probability that p is accepted is at least 

(1 - e - ^=)(1 - d^poly(L)) > 1 - 6 - -L _ d^poly(L) 
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Now we have to show that 

-L + d^polyiL) 

can be made smaller than 5. Choose k such that 1/yk is smaller than 5/2, i.e., 
k > 16/<5 4 . 

As already noted, it is possible to choose d such that (f 1 / 4 poly(L) < 5/2. To see 
this, recall that L is here not the length of the original circuit but of the meta circuit 
U. However, it increases only polynomially in k and the least possible value of a is 
0(2 -v/ ^/ ln2 ) decreasing faster than the reciprocal of every polynomial. Therefore d 
can be chosen as required. □ 



4 Which temperature is required 
for solving FQMA and QMA? 

Obviously low energy state preparation can be achieved by cooling. But clearly the 
invention of efficient cooling mechanism is a highly non-trivial problem for theoreti- 
cians and experimentalists. In the context of quantum computation algorithms have 
been suggested to prepare low energy states on a system with universal quantum com- 
putation capabilities [2J El • At first sight there seems to be a fundamental difference 
between the task to cool a system with given Hamiltonian and to prepare a thermal 
equilibrium state according to a Hamiltonian that is specified by classical information 
only. In the latter case the Hamiltonian is only virtual, i.e., not physically present. 
However, from the point of view of complexity theory this does not make any differ- 
ence: It seems that every cooling mechanism can be translated into an algorithm that 
prepares low temperature states on a quantum computer: Given the postulate that 
every physical process can efficiently be simulated on a quantum computer one can 
clearly simulate the time evolution according to the Hamiltonian of the system that is 
to be cooled and the interaction between this system and its environment (the entropy 
sink). Whether each algorithm can be translated into a cooling mechanism depends 
on the control operations which are available. Consider for instance an n + /c-qubit 
system with pair-interactions forming a connected graph. Assume furthermore that 
all one-qubit transformations can be implemented on a time scale that is smaller than 
the time evolution according to the interaction. Then decoupling-techniques can be 
used to "switch off" unwanted interactions or transform them to different interaction 
types in such a way that universal quantum computation is possible (e.g. [T3 l IT U ITB] ) . 
Hence one can implement every quantum algorithm which transports entropy from the 
n qubits to the k ancilla qubits. 

The efficiency of algorithms preparing equilibrium states seems to depend strongly 
on the demanded temperature |12[ |2j and also the required thermodynamical resources 
increase strongly with temperature tending to the absolute zero |15j . 

We consider the question which temperature is needed to solve FQMA. Let us first 
consider Low Energy State Preparation (see Definition 13}. The gap between the energy 
value d n to be achieved and the ground state energy g n is essential for the required 
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temperature. To see this define a cutoff energy value by 

1 

C n . — CL n -\- 2' n Q"n) • 

In the following we shall drop the index n. Consider the Gibbs equilibrium state p with 
absolute temperature T. Let be a basis of eigenstates of H n with corresponding 

eigenvalues (Ej)j. Then we have 

3 

with 

e -E s /{KT) 
P 3 := £ e -E 3 /{KT) ' 

where K is the Boltzmann constant. We estimate the probability pj for all j with 
Ej > c. We have 

e -E 3 /{KT) 1 
Pj < e -9/(KT) + e -Ej/{KT) ~ e (c-a)/(KT) 

We have at most 2 n — 1 < 2 n states above c. Therefore we can estimate the mean 
energy of p as follows: 

tr(pH) = ^P 3 Ej<c + 2 n e(c J /(KT) E max 
3 

1 gn In 2 

= a + -{d-a) + e{d _ a)/(2KT) Emax , (3) 

where E max is the largest eigenvalue of i7. It is 0(n 3 ) since the Hamiltonian contains 
n(n — l)(n — 2)/6 three-body interaction terms. The mean energy is asymptotically 
less than d n if (d — a)/(2KT) > nln2. This can be achieved by the temperature 

Tn < Kn2ln2q(n) ~' Kq(n) ' ^ 

i.e., the temperature must only decrease as the reciprocal of a polynomial in n. 

If q is chosen appropriately, the temperature in eq. (J5J) is sufficient to solve FQMA- 
problems provided that n is the size of the register where the meta circuit U acts on 
plus the size of clock register. 

Note that a cooling mechanism is in some sense less worth than an extremely good 
refrigerator. The cooling mechanism required to solve FQMA has only to be able to 
bring the system in an extremely low temperature after one has specified the Hamil- 
tonian. It is comparable to a fridge with some buttons to provide it with information 
whether cheese, SAT or ice cream should be cooled, i.e., one specifies 



Hcheese Or H SAT OT H. 



ice cream 



A more explicit upper bound on the required temperature can easily be derived for 
solving QMA-problems instead of FQMA. Assume that the circuit either accepts all 
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states with probability less than e or there is a state that is accepted with probability 
1 — e. Let L be the length of the circuit and H be the Hamiltonian corresponding 
to U. We assume that all either there is an eigenvalue of smaller than e/(L + 1) 
or that all eigenvalues are larger than (1 — e)/(L + 1). This can always be achieved 
by amplification and considering instead the Hamiltonian corresponding to the meta 
circuit. The number of qubits the Hamiltonian acts on is given by n := n x + m x + L + 1 
where U acts on an n x + m x dimensional register as in Definition^ We choose T n with 

T l ~ 2€ fK\ 

ln< 4ln2K(L + l)n- W 



Set the cutoff value at 



l + 2e 



4(L + 1) 

Then one can check using eq. @ that if there is a witness then for large n the mean 
energy is smaller than 

d= 2(1+1) " 

We conclude that for T n satisfying inequality © either the mean energy is above 
(1 — e)/(L + 1) or below 1/(2(L + 1) depending on the existence of a witness. Since 
the gap between these values is greater than the reciprocal of a polynomial they can 
efficiently be distinguished using some copies of the same Gibbs state. 

Let us just mention a similar bound on the temperature which is required to find 
the solution of the classical (NP-complete) Independent Set Problem. In the Ising- 
Hamiltonian on n spins considered in ED] the energy value of the ground states encodes 
a graph theoretical NP-problem. The energy gap between the ground state and the 
first excited state is a constant value AE. Then 

T< 1 



In 2KAEn 



is sufficient to have for large n an almost sure answer. 
This can be seen by choosing the cutoff value 

c:= a + -AE . 

4 

If the mean energy of the system is above 

1 * ^ 
a + -AE 
2 

we can conclude that the answer of the NP-problem is "no" , otherwise it is "yes" . 
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